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ABSTRACT 

Earlier studies of grain alignment dealt mostly with interstellar grains that have strong internal 
relaxation of energy which aligns grain axis of maximum moment of inertia with respect to rain's 
angular momentum. In this paper, we study the alignment by radiative torques for large irregular 
grains, e.g., grains in accretion disks, for which internal relaxation is subdominant. We use both 
numerical calculations and the analytical model of a helical grain introduced by us earlier. We 
demonstrate that grains in such a regime exhibit more complex dynamics. In particular, if initially the 
grain axis of maximum moment of inertia makes a small angle with angular momentum, then radiative 
torques can align the grain axis of maximum moment of inertia with angular momentum, and both 
axis of maximum moment of inertia and angular momentum are aligned with the magnetic field when 
attractors with high angular momentum (high-J attractors) are available. For the alignment without 
high- J attractors, beside the earlier studied attractors with low angular momentum (low- J attractors), 
there appears new low-J attractors. The former and later cases correspond to the alignment with long 
axes perpendicular and parallel to the angular momentum, respectively. In addition, we also study 
the alignment of grains in the presence of strong internal relaxation, but induced not by a radiation 
beam as in earlier studies, instead, induced by a complex radiation field, that can be decomposed into 
dipole and quadrupole components. We found that in this situation, the parameter space q rnax : for 
the existence of high-J attractors in trajectory maps is more extended, which entails higher degrees of 
polarization expected. Our obtained results are useful for modeling polarization arising from aligned 
grains in molecular clouds and accretion disks. 

Subject headings: Polarization -dust extinction -ISM: magnetic fields 



1. INTRODUCTION 

Magnetic fields play a crucial role in many astrophys- 
ical processes (e.g., star formation, accretion disks, cos- 
mic ray transport). An important and easily available 
way to map magnetic fields is to observe the polarization 
of emission or absorption by aligned grains (Goodman 
et al. 1995; Hildebrand et al. 2000, 2002; Crutcher et 
al. 2004). The reliable interpretation of polarimetry in 
terms of magnetic fields requires, however, a solid un- 
derstanding of the alignment of grains. This is the ma- 
jor motivation of work in the field of grain alignment. 
An additional motivation to understand when grains get 
aligned comes from the necessity of separating the po- 
larization by dust grains from the polarized Cosmic Mi- 
crowave Background (CMB) signal (see Lazarian 2008). 

Grain alignment is a complex area of studies related to 
many physical effects of different time scales (see review 
by Lazarian 2007). It is also the research area associated 
with the names of Lyman Spitzer and Ed Purcell, who 
decisively contributed to understanding of a number of 
fundamental physical processes of grain alignment. 

A number of mechanism has been known to induce 
grain alignment (see Lazarian 2007). However, recently 
grain helicity has been discussed as the necessary element 
of successful alignment. Note that the idea of grains 
having twist and therefore interacting differently with 
left and right polarized photons can be traced back to the 
pioneering work by Dolginov & Mytrophanov (1976). A 
more recent work in the area showed that a similar effect 
is present when irregular grains interact with a gaseous 



flow (Lazarian & Hoang 2007b). 1 

In this paper we concentrate on studying grain align- 
ment by radiative torques (hereafter RATs). However, 
the intrinsic similarity of the radiative and mechanical 
torques acting on helical grains (see discussion in Lazar- 
ian & Hoang 2007b) makes our present results applicable 
to the alignment by mechanical torques. Note that for 
the alignment of helical grains, the magnetic field acts 
only through inducing the Larmor precession. Therefore 
the alignment may potentially occur with both longer 
and shorter grain axes perpendicular to the magnetic 
field. Observations of interstellar dust indicate that the 
alignment with the long axes perpendicular to the mag- 
netic field is dominant in the diffuse gas. Therefore we 
shall call this alignment right alignment, and the opposite 
type of alignment we shall call wrong alignment. Differ- 
ent alignment processes produce either right or wrong 
alignment and the goal of the grain alignment theory is 
to predict correctly both the direction and the degree of 
the expected alignment. 

As dust grain scatters or absorbs photons it experi- 
ences radiative torques. These torques can be stochastic 
or regular. Stochastic RATs arise from, for instance, a 
spheroidal grain randomly emitting or absorbing pho- 
tons. The latter process, for instance, was invoked by 
Harwit (1970) in his model of grain alignment based 
on grains being preferentially spun up in the direction 
perpendicular to the photon beam. However, Purcell 

1 Such flows arise naturally in the reference of a grain as the 
grain interacts with ubiquitous MHD turbulence (Lazarian & Yan 
2002; Yan & Lazarian 2003; Yan, Lazarian, & Draine 2004). 
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& Spitzer (1971) showed that the randomization arising 
from the same grain emitting thermal photons makes the 
achievable degree of grain alignment negligible. In fact, 
they showed that stochastic RATs arising from grain 
thermal emission is an important process of grain ran- 
domization. More recently, grain thermal emission was 
analyzed as a source of the excitation of grain rotation as 
well as its damping in relation to the rotation of tiny spin- 
ning grains that are likely to be responsible for the so- 
called anomalous foreground emission (Draine & Lazar- 
ian 1998). 

RATs were first discussed by Dolginov (1972) in terms 
of chiral, e.g. hypothetical quartz grains. The idea is 
that starlight passing through such grains would spin 
them up. Later, Dolginov & Mytrophanov (1976) con- 
sidered an irregular grain model of two spheroids twisted 
to each other, made of more accepted materials, e.g. sil- 
icate grains, and claimed that these grains will be both 
spun up and aligned by regular RATs. This work was, 
unfortunately, mostly ignored for 20 years. 

Grain alignment by RATs drew much attention when 
Bruce Draine modified the Discrete Dipole Approxima- 
tion Code(Draine & Flatau 2004, hereafter DDSCAT) 
and enabled later researchers to calculate RATs for ir- 
regular grains. 2 

Draine & Weingartner (1996, hereafter DW96) treated 
RATs kind of pinwheel torque, and they found 

that RATs can accelerate grains to suprathermal rota- 
tions. Later, Draine & Weingartner (1997, hereafter 
DW97) introduced many essential elements of the mod- 
ern treatment of RATs, for instance, phase trajectory 
maps. They confirmed Dolginov & Mytraphanov (1976) 
hypothesis that RATs can produce the alignment by their 
own. However, the important moments of grain dynam- 
ics, so-called crossovers, were not treated correctly in 
the paper. Crossovers, are special events in the dynam- 
ics of grains rotating subject to regular torques. Dur- 
ing crossovers grain angular velocity approaches its min- 
imal value and models of crossovers caused by pinwheel 
torques in Spitzer & McGlynn (1979) and in Lazarian 
& Draine (1997) envisaged grain flipping with the direc- 
tion of the angular momentum staying the same. DW97, 
however, assumed that the RAT crossovers are different 
and during the RAT-induced crossovers angular momen- 
tum changes its direction to the opposite. This treatment 
of crossovers resulted in grains having cyclic trajectories, 
which were, in fact, an artifact of the adopted treatment, 
as was shown later (Lazarian & Hoang 2007a, henceforth 
LH07a). In fact, LH07a showed that instead of cyclic 
trajectories, one gets situations when RATs tend to slow 
grains down. 

A more general analysis for dynamics of RAT align- 
ment was done in Weingartner & Draine (2003, hence- 
forth WD03), where thermal fluctuations within grain 
were accounted for. These fluctuations were shown by 
Lazarian (1994, see also Lazarian & Roberge 1997) to 
induce grain thermal wobbling, which amplitude gets 
larger as the grain angular momentum approaches its 
thermal values. For their choice of parameters, WD03 
found that an attractor with low angular momentum, so- 

2 The only work on RATs before that was Lazarian (1995) which 
corrected some of the points in the Dolginov & Mytrophanov (1976) 
study, but underestimated the importance of RATs. 



called, low-J attractor, appears on the trajectory phase 
maps. 3 The study, however, was limited in terms of pa- 
rameter space explored. In fact, it was a numerical study 
of alignment of a single grain subjected to the electro- 
magnetic wave of a single frequency at a single incident 
direction. Therefore implications of the study for the 
problem of the RAT alignment were difficult to evaluate. 

We feel that the numerical work above done with 
DDSCAT should be treated as an experimental attack on 
the problem of the RAT alignment. The limitations of 
such an approach are self-evident. For instance, to obtain 
predictions one should analyze many hundreds of differ- 
ent grain shapes and many wavelengths and many inci- 
dent directions. The study in DW97 included 3 shapes 
and resulting torques looked very different, which made 
one wonder what causes the alignment 4 . At the same 
time, unfortunately, the analytical study in Dolginov & 
Mytrophanov (1976) was in error. 

The first analytical model (henceforth AMO) of RATs 
consistent with the numerical calculations was proposed 
in LH07a. This study provided simple analytical expres- 
sions for the RATs (see also §2). In particular, the func- 
tional forms of the torques acting on grains of differ- 
ent shapes are similar, with the difference between the 
grains of different shapes amounting to a single parame- 
ter, termed q max -r&tio, which is the ratio of the magni- 
tudes of two first components of torques in the grain sym- 
metry system of reference where the radiation direction is 
the symmetry axis. AMO allowed theoretical predictions 
and the analytical treatment of the RAT alignment. As 
a result, studies of a parameter space seized to be an in- 
surmountable problem. In LH07a, similar to DW97, the 
thermal fluctuations were intentionally disregarded and 
in such a setup the low-J attractors, rather than cyclic 
trajectories, were reported. In fact, it became clear that 
RATs, over a large part of the parameter space, spin 
down rather than spin up grains. Using AMO, LH07a 
obtained the values q max for which grains have only low- 
J attractors and when they have both low-J and high- J 
attractors simultaneously. 

AMO was extensively used in all our papers that fol- 
lowed. For instance, in Hoang & Lazarian (2008a, hence- 
forth HL08a), we studied the RAT alignment in the 
presence of thermal fluctuations using both AMO and 
DDSCAT for an extended sample of grain shapes, radi- 
ation direction and wavelength. We found that irregu- 
lar grains do not stop completely, but rotate at a rate, 
which is comparable with the rate of thermal rotation 
at the dust temperature, which agrees with an earlier 
example in WD03. Importantly, HL08a considered the 
RAT alignment in the presence of gaseous bombardment 
and reported a new effect, namely, the transfer between 
the low-J and high- J attractors, in situations when the 
both low-J and high- J attractors were present simulta- 
neously 5 . Thus, counterintuitively, the collisions were 

3 Conventionally, an attractor with angular momentum J larger 
than the thermal value is called a high- J attractor, and the attrac- 
tor with J of the order of thermal value is called low-J attractor. 

In relation to this one of us (AL) recalls that Lyman Spitzer 
after studying the work on RATs suggested that one should seek 
for some simple trigonometric fits to torques that would produce 
the alignment. 

5 Note that the parameter space for the existence of low-J and 
high- J attractors are similar and known through the LH07a study. 
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found to increase the alignment. 

The most important practical implication of the LH07a 
and HL08a studies was that, over a large range of the 
parameter space of g™ la:l: -ratio and angles between the 
magnetic field and radiation direction, the grains rotate 
thermally and therefore, as a result the aforementioned 
thermal fluctuations of grains, and exhibit the reduced 
alignment of grain axes with respect to J. This opened 
possibilities of estimating the expected alignment and 
comparing it to the alignment inferred from observations. 
Importantly enough, Lazarian & Hoang (2008, hence- 
forth LH08) reported, that in the presence of superpara- 
magnetic inclusions (Jones & Spitzer 1967, Mathis 1986, 
Bradley 1994, Martin 1994, Goodman & Whittet 1994), 
the high- J attractors always exist. This finding entails 
a very non-trivial effect, namely, that, in the presence 
of both superparamagnetic inclusions, grains, that oth- 
erwise would rotate subthermally, get into the state of 
fast suprathermal rotation. According to HL08a, this 
also means that all grains eventually get into the state 
of fast rotation and perfect alignment. From the point of 
view of observations, this opens prospects of testing the 
existence of superparamagnetic inclusions by measuring 
the degree of polarization. 

This paper continues our studies in LH07a and HL08a. 
The major issues we address below are (1) the effect of 
internal relaxation, and (2) the effect of complex radia- 
tion field on the grain alignment. 

The first issue is important in terms of the alignment 
of large grains, e.g. grains within proto-planetary accre- 
tion disks. The alignment of such grains was assumed in 
a recent modeling of polarization from T-Tauri accretion 
disks in Cho & Lazarian (2007). Grains which are im- 
portant for such a modeling may, however, should be suf- 
ficiently large that the internal relaxation time for them 
is longer than the alignment time. Therefore, how grains 
get aligned in the absence of efficient internal relaxation 
is an important question. 

The second issue is vital for modeling the polarization 
from both disks and dark clouds. The latter modeling 
which used the RAT alignment was undertaken in Cho 
& Lazarian (2005), Pelkonen et al. (2007) and Bethell 
et al. (2007). All these studies assumed rather naive 
models of alignment. In particular, the actual modeling 
would require taking into account the actual complex 
radiation field that is being experienced by the grain, 
rather than to assume that the radiation is coming from 
a point source. 

In what follows, in §2 we briefly describe the AMO, 
grain torque- free motion, and calculations of RATs. In 
§3 we compare the Barnett and nuclear relaxation times 
with the radiative alignment time and identify grain size 
when the internal relaxation can be disregarded. In §4 
we study the properties of RATs and the resulting RAT 
alignment for large grains in the absence of the internal 
relaxation. Grain alignment by RATs arising from dipolc 
and quadrupole components of radiation fields is studied 
in §5. We discuss our findings in §6, and our summary 
is presented in §7. 

2. RADIATIVE TORQUES 

2.1. Definition of radiative torques 

Considering only the anisotropic component of the ra- 
diation field, the radiative torque resulting from the in- 




FlG. 1. — Upper panel: Orientation of a grain, described by three 
principal axes Oi, 0,2 and S3, in the lab coordinate system §1, §2, e~3 
is defined by three angles O, (3 and The direction of incident 
radiation beam k is along ei. Lower panel: The orientation of 
angular momentum J in the lab coordinate system, tp is the angle 
between the magnetic field B and k, £ is the angle between J and 
B, and <f> is the Larmor precession angle of J about B. 

teraction of a radiation beam with a grain is then defined 

by 

r rad = ^^ 7 q r (e, /?,$), (i) 

where a is the grain size, 7 is the degree of anisotropy, A 
and u ra( i are the mean wave length and total energy den- 
sity of the radiation field, respectively; Qr is the RAT 
efficiency vector depending on angles 0, (3 and <&, and its 
overline denotes averaging over the spectrum of radiation 
field (DW96); 8 is the angle between the axis of max- 
imum moment of inertia ai and the radiation direction 
k, (3 is the rotation angle about ai and $ is the preces- 
sion angle of ai about k (see Fig. \\^ipper). In the lab 
coordinate system eie^e^, Qp is decomposed as 

Qr(6, /3, *) = Q ei (9, (3, *)ei + Q e2 (6, {3, $)e 2 

+Q e3 (e,/?,$)e 3 . (2) 
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Fig. 2. — Sketch of an AMO consisting of a weightless mirror 
attached to a perfectly reflecting ellipsoid. Principal axes Si , (12 
and (13 define a body coordinate system, and ei,e2 and £3 define 
the lab coordinate system. is the angle between a\ and the 
radiation direction k||ei . N is the normal vector lying in the plane 
0,20,3, and a is the tilted angle of the mirror in the body coordinate 
system, r is the radius vector of length Zi connecting the ellipsoid 
to the mirror, I2 is the mirror size and l\ 3> I2 is assumed. 

2.2. Analytical model of RATs 

LH07a proposed an analytical model of RATs arising 
from the interaction of a helical grain, consisting of a 
mirror attached to an ellipsoid body with a radiation 
beam. Although this model was derived in the geometric 
optics approximation, LH07a proved that it represents 
well RATs obtained by DDSCAT. 

The geometry of AMO is depicted in Figure [2] It has a 
weightless mirror tilted by an angle a about the axis &2 
attached to an ellipsoid. Detail comparisons of RATs on 
the model and irregular grains suggested a = 7r/4. For a 
given direction of the grain with respect to the radiation 
direction k, RATs arising from the scattering of photons 
on the mirror are functions of the angles G,/3 and $. 
LH07a calculated three components of RATs at $ = 0: 
g ei (e,/3,0),Q e2 (e,/3,0) and Q es (e,/3,0). 

LH07a found that whether the RAT alignment has 
high-J or low-J attractors depends on the ratio of torque 
components q max = Q™ l ® x / Q™ ax (see Appendix A). 
Therefore, this parameter has been used to character- 
ize the RAT alignment in LH07a. For irregular grains, 
LH07a found that q max depends on the wavelength of ra- 
diation, grain shape, size and composition, ranging from 
to 10 2 . Thus, we treat q max as a variable parameter 
in this paper. 

To conform with numerical calculations using 
DDSCAT, we combine the functional forms of RATs 
from AMO with the scaling for the magnitude of RATs 
as a function of wavelength and grain size. The final 
RATs used in our series of papers are given in Appendix 
A. We note that for the AMO, the functional forms of 
RATs do not depend on wavelength, so that we can 
simplify the notations by using Qp instead of Q r . 

2.3. Torque-free motion 

The motion of a triaxial grain in the absence of external 
torques, namely, torque-free motion, is well-known (see 
Landau & Lifshitz 1976; also WD03). Let consider an 
irregular grain with angular momentum </, and rotational 




Fig. 3. — Angular velocity components as functions of time t/Pt 
with Pt being the rotation period, for different values of p for an 
irregular grain with I\ : 1% : I3 = 1 : 1.5 : 1 and in a positive 
flipping state, i.e., ui\ and L03 both being positive for p < I\/l2- 
The value of p corresponding to each curve is labelled. 

energy E. It is convenient to introduce a dimensionless 
parameter 6 



P 



2hE 



(3) 



where I± is the inertia moment about the axis of maxi- 
mum moment of inertia ai, and other inertia moments 
are smaller than I±, i.e., I± > I2 > la- Here p spans in 
the range 1 to -Z1//3, where the lower and upper limits 
correspond to the grain rotation about the axis of maxi- 
mum moment of inertia ai and axis of minimum moment 
of inertia a3, respectively. 

For the torque-free motion, the angular velocity com- 
ponents about the three principal axes, u>i , 0J2 and u>3 can 
be obtained analytically as functions of p and time t (see 
WD03; HL08a for their expressions). In Figure [3] we plot 
wi, 0J2 and LU3 as functions of time t/ Pt with P t being the 
rotation period, for different p for an irregular grain with 
ratio of inertia moments I\ : 1% : ^3 = 2 : 1.5 : 1. The 
plot corresponds to a positive flipping state, i.e., u>\ and 
LU3 are both positive for p < I\/l2- Similarly, when lo\ 
and W3 are both negative for p < I1/I2, we call negative 
flipping state (see HL08a for details). 

We can see that when p — > 1 or I1/I3 = 2, the rotation 
is stable about ai and a3 axes, respectively. However, 
the rotation is not stable about the intermediate axis a2 
when p = 1.32 and 1.34, which are very close to I1/I2 — 
1.333. In fact, the middle panel shows that as p — > I1/I2, 
the grain gets into a negative state with I20J2/ J = — 1 for 
almost one half period, and then reverses its rotational 

6 This parameter was denoted by q in WD03, however, we use p 
to avoid the possible confusion with q max . 
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state, staying there for another half period. Of course, 
at p = I1/I2, we have I2OJ2/ J = 1 and u)\ = UJ3 = 0. 

Because the timescale of torque-free motion is much 
shorter than the internal relaxation timescale and other 
dynamical timescales (e.g., Larmor precession, gas damp- 
ing times), we can average RATs over such a motion 
to obtain Qr(9, (3, 3?) where the tilde denotes averaging 
over torque-free motion. Resulting RATs are functions 
of £, 4>, ip,p and J. 

2.4. Radiative torques in alignment coordinate system 

To study the alignment of the angular momentum J 
with the magnetic field B, we represent RATs in the 
spherical coordinate system J,£,<f>, the so-called align- 
ment coordinate system (see the lower panel in Fig. [I]). 

In this alignment coordinate system, the RAT T ra d can 
be written as 

Trad = M[F(t, </>, p, J)i + G(£, 0, p, J)$ + Hfo 0, p, J)&) 

where the dependence on the angle between the radiation 
direction k and B, ip, is omitted, 



F, which is the torque component parallel to £, acts to 
change the orientation of J with respect to B. H, the 
component parallel to J, is to spin grains up and G in- 
duces the precession of J about the magnetic field or 
radiation. The later are given by 

(p,p, J) =Q ei (£> <^)P) J) (cos -0cos £ — sin ^sin £cos (p) 

+0e 2 (£> 4»:P> ./)(sim/)Cos £ + cos ^sin £cos <j>) 

+Qe 3 (£, <P, P, J)sin £sin (p, (6) 

F(£, (j),p, J) —Q ei {£,, 4>,P, J) (—sin -0cos £cos <f> — sin £cos ip) 

+Qe 2 (£j §-,Pi J) (cos V'cos £cos <p — sin £sin ip) 

+Qe 3 (£, <P, P, ^)cos £sin <f>, (7) 

where £, ^ and (p are angles describing the direction of J 
in the lab coordinate system (see Fig. [TJ lower panel). 
Here Q ei (£, 0,p, J) with i — 1, 2 and 3 are RAT compo- 
nents along ei axes, respectively, obtained by averaging 
Q ei over torque-free motion. 

We define also the torque components along grain axes 
ai , a.2 and a3 as 

Q Oi (£,0,p,J) = Qr-^for i=l to 3, (8) 

and thus, the work done by RATs per second 

Qu, = Qr-^ = Q ai 0Jl + Qa 2 U2 + Qa 3 UJ 3 . (9) 

Here u>i for i = 1, 2 and 3 corresponds to the components 
of angular velocity along three grain axes a^ Their ex- 
pressions for torque-free motion are given in Appendix 
A of HL08a. In addition, from Figure [TJ we have 

ai = eiCos9 + e2sin9cos$ + e3sin8sin(f>, (10) 
a 2 = — eisin9cos/3 + e 2 (cos9cos$cos/3 — sin<I>sin/3) 

+e 3 (cos9sin$cos/3 + cos$sin/3) , (11) 
a3 =eisin9sin/3 — e 2 (cos9cos$sin/3 + sin<E>cos/3) 

+e 3 (— cos0sin$sin/3 + cos$cos/3). (12) 




a 3 

Fig. 4. — Sketch of the irregular grain, shape 1, built from 
13 cubes adopted in our calculations. The inertia moments along 
principal axes are I\ : Ii : I3 = 1.745 : 1.61 : 0.876. 



2.5. Calculations of RATs 

Taking analytical expressions Q ei (Q, [3,$ = 0) with 
i = 1,2 and 3 from AMO (see Appendix A), we first cal- 
culate Q ei (Q,(3,$) by using equations (|A13[> - (|A15|) for 
9 S [0,7r], /3 S [0, 2?r], and $ G [0, 2?r]. Next, we calculate 
RATs for J from 0.1 to 2QI x u) T , £ from to tt, <p from to 
27r, and p from 1 to I\ /I3 and for a given light direction ip. 
For a given set of parameters J, £, <p, p, we use coordinate 
transformations to get angles 9, f3 and $ as functions of 
Euler angles (see Appendix A in HL08a). Finally, torque 
components Q ei (£, ip, (p) are obtained using interpolation 
for Q ei (9, (3, <l>). Then we average Q ei (£, ip, <p) over Eu- 
ler angles, i.e., averaging over torque- free motion, to get 

Substituting Q ei (^,ip,(p)} into equations ((6]), ([7]) and 
© we obtain F(£,ip,(p),H{€,ip,<p) and Q^,ip,4>). We 
then average them over angle to get the Larmor preces- 
sion averaged values. Since these torques are functions of 
J, £,p, we denote them by (F) <j,(J,£,p), (H)$(J,£,p) and 
(Qu>)cj>(J-)£,p)-) where the dependence on ip is omitted. 
We will use these results to solve equations of motion for 
the alignment in §4. 

3. TIME SCALES: GAS DAMPING, INTERNAL 
RELAXATION AND RADIATIVE ALIGNMENT 

3.1. Gas damping time 

We represent the magnitude of torques in terms of the 
thermal angular momentum, and the gas damping time. 
The later, for the sake of simplicity, can be obtained for 
an oblate spheroid, with the moments of inertia Iu = 
Ii, 7j_ = I2 = I3J The expressions for triaxial ellipsoids 

7 For the oblate spheroid with major and minor axes a and b, 
7|| = 2Ma 2 /5 = (87rp/15)6a 4 , and I± = (4vrp/15)a 2 6(a 2 + b 2 ), but 
sometimes we define the equivalent sphere with the same volume 
of the grain, and use Ii = (8ir/ 45) po^a 5 with ai for i = 1,2 and 3 
being dimensionless factors (see WD03). 
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are mor complex, but do not change results substantially. 
Then, the thermal angular momentum is given by 



Jth — \/ I\\kBT„a S — 



87rpa 5 s 
15 



kBTna 



gas • 



5.89 x 10 



-20„5/_2 -1/2 -1/2 j,l/2 



pas " 



(13) 



with p density of material within the grain and p 



jo/3 gem , a grain size, a_5 



a/10" 



cm, s 



s/0.5 



b/a and T gas gas temperature and T gas 



with i 

2gas/100K. And the damping time due to gas colli- 
sion (see Roberge, DeGraff & Flatherty 1993, hereafter 
RDF93) is 

t 3 J|1 

°gas — 



- 2.3 x 10 12 a_ 5 sf - 1 / 2 



40r n H m H a' l v th r\\ 
30 cm -3. 



(14) 



where v t h is the thermal velocity of a gas particle with 
density uh, rnu is the hydrogen mass, Tu is the geomet- 
rical parameter which is unity for sphere (see RDF93), 
and we adopted standard parameters in the ISM (see 
Table 1 in HL08a). 



3.3. Radiative alignment timescale 

The characteristic timescale for RATs to accelerate the 
grain from J = to J — I\to is estimated by 

I\oj 1\uj 



trad 



rad 



M(H) 



(18) 



where (H)j, is the spin-up component of RATs, that is 
averaged over the torque-free motion and the precession 
angle 4> using equation ([5]). 

Using M given by equation ([5]) and I\ = 16pa 5 s/3 for 
the brick of 2a x 2a x 2b sides, we obtain 



(19) 



t _ 2 3 / 2 (16pfcr gas /3) 1 / 2 a 1 /V/ 2 
Aw rad 7(ff)0 

For standard parameters of the ISM, it yields 

+ ^ O0v in io-T " 1/2 1-2 p,m uisrf 10~ 3 , . 
t ra d ~ 2.8 x 10 'pT gas saJ_ b = /m ^0) 



A U r ad l{H) 



<t> 



where s — s/0.5 with s = b/a being the ratio of long 
to short axes, uisrf is the energy density of interstellar 
radiation field (hereafter ISRF; see Mathis, Merger, & 
Panagia 1983). 



3.2. Internal relaxation timescales 

A rotating paramagnetic grain experiences internal dis- 
sipation of energy due to Barnett and nuclear relaxation 
(see Purcell 1979 and Lazarian & Drain 1999b, hereafter 
LD99b, respectively). In addition, superparamagnetic 
inclusions can increase significantly the rate of internal 
relaxation (Lazarian & Hoang 2008). 

Disregarding other relaxation processes within the 
grain, e.g., inelastic relaxation (Purcell 1979; Lazarian 
& Efroimsky 1998) as well as the relaxation due to su- 
perparamagnetic inclusions, the internal relaxation rate 
arising from Barnett and nuclear relaxation for a brick 
of 2a x 2a x 26 sides (see Fig. [TSJ lower panel) is then 



tl 



L Bar ^ h nucV 



: 1.310- 7 /r 2 s[0.5 + O.UhffaZlsfgasf^ 1 



J 



(h-l)f(J). 



(15) 



J\U>T , 

Here the Barnett and nuclear relaxation times tsar and 
tnuci are taken from Hoang & Lazarian (2008b; HL08b), 

1 /2 

s = b/a, h = 7||/Jj_, lot = {^kT gas / I\j is the thermal 
angular velocity of grain at the gas temperature T gas , 

Td = 7rf/15K with Td being the dust temperature, and 



/(J) = 1.610 5 [1 



L0\T n 



0.47[1 



2 J 1 V 2 

For grains larger than 1/mi, f( J) w 1.6 x 10 5 for J = J t h- 
In above equations, io\ — Jcos9/Iu, and r n is given by 



T n — l/( T ne + T nn) 



(17) 



where r ne = 3 x 10" 4 (2.7 /g n f (l0 22 cm- 3 /n e ) s and 
T nn = «/(3.8g„n„p„) w 0.58T„e (n e /n n ) with electron 
and nucleus density n e and n n are the relaxation time 
of interaction nucleus-electron and nucleus-nucleus spins, 
respectively. Also, r e ; = 2.9 x 10~ n s is the spin-spin re- 
laxation time of electronic spins, and cos# = 1/2 is chosen 
(see LD99b; HL08b) 



3.4. Radiative alignment time versus internal relaxation 

time 

To compare the effect of RATs and internal relaxation, 
let us estimate the ratio of their timescales. Following 
equations (TT5")) and (|2"0|) we obtain for a brick grain 



3.7 x 10 3 p- 15 [0.5 + 0.125s 2 ] 5 



T, 



1.5 



1 



10 



-3 



- ) f(J), 



Td U ra dX 

(21) 



Here T d = r rf /15K, u rad = u rad /u IS RF, and A = 
A/1.2 pm. uisrf = 8.64 x 10 _13 erg cm -3 is the energy 
density of interstellar radiation field (hereafter ISRF; see 
Mathis, Mezger & Panagia 1983), and 7 = 0.1. 
Using the scaling of RAT magnitude obtained in 



LH07a: (77 

get 

^rad 



\Qr\ 



0.4 



for A > 1.8a. We 



Thl A 2 



8.4 x 10 5 ^- 1 - 5 [0.5 + 0.125S 2 ] 2 aI§- 8 gas 



I — I f(J)- 



For A < 1. 



(77) 



0.4, 



trad - 

- — = 70p 



JTl.S 1 

L5 jq~ 6,5 - 



(22) 

— ) f{jy (23) 



Both the increase of the grain size and the radiation in- 
tensity lead to the decrease of t ra d- 

For a grain size a = 1.2 pm, f(J = J t h) ~ 10 5 , so it 
can been seen from equation (|23p that t ra d/ti n t ~ 0.67 
for the ISRF (see also Fig. HJ. 

Figure [5] shows t gas , ti nt and t ra d as functions of a for 
the ISM obtained from equations Q3}, (TS]) and ([2Ti]). re- 
spectively. It can be seen that for grains smaller than 
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Fig. 5. — Timescales as functions of grain size for internal relax- 
ation (ti nt ), radiative (t rat j) and gas damping (t sas ) calculated for 
J = J t h and the ISM. Filled circles denote where t; nt = t rac j. 

~ 1/im, ti n t/t rac i <C 1, and decreases steeply with a de- 
creasing. For grains larger than ~ 1.2pm, we see that 
Unt /Uad > 1, and tint /t rac i increases rapidly with a. 
When the mean energy density u ra d increases bylO 2 time, 
the size corresponding to tint = t ra d decreases from 1.5 
to 0.6/xm. 

Note that earlier works on grain alignment dealt with 
the alignment of interstellar grains with size in the range 
from 0.005 to 0,25pm. For this range of grain size, in- 
ternal relaxation is very strong, so that the average of 
RATs over thermal fluctuations arising from internal re- 
laxation was accounted for (WD03; HL08a). In some 
circumstances, e.g. accretion disks, molecular clouds, 
where larger grains corresponding to weak internal re- 
laxation are expected, we need to study the internal and 
external alignment at the same time. 

Another characteristic timescale involving in the grain 
dynamics is the Larmor precession time of grain magnetic 
moment about an ambient magnetic field. Due to the 
Barnett effect, a rotating paramagnetic grain develops 
a magnetic moment psar which is proportional to the 
angular velocity (see Dolginov & Mytraphanov 1976). 
The value of pBar is given by 



_ x(Q)vn 

PBar — U), 

9VB 



(24) 



where x(0) = 4.2 x Wr 2 f p T^ x with f p being the fraction 
of paramagnetic material, is the magnetic susceptibility 
at zero frequency, V is the volume, g is gyromagnetic 
ration, which is ~ 2 for electrons, and \ib — e/2m e is 
the Bohr magneton. 

The Larmor precession time in an external magnetic 
field B is then 



tr, = 



2ttIiuj 



A*Bc 



B 



4.2 x 10V, 



pT d 
'BX(O) 



(25) 



where B = B/(5pG) and x(0) = x(0)/10 -3 , and we used 
h = 8wpa 5 /15 for inertia moment. 

It is easy to see that this time scale is shorter than both 
internal relaxation and radiative alignment timescales for 
grains larger than about /jm (see Fig. [5]), so that we can 
average RATs over the Larmor precession while dealing 
with the overall alignment. 



4. RAT ALIGNMENT IN THE ABSENCE OF INTERNAL 
RELAXATION 

Now let us study grain dynamics for grains larger than 
lpm in the diffuse interstellar medium. For this range 
of grain size, U n t > trad, the influence of internal relax- 
ation to grain alignment can be disregarded. As a result, 
we follow both the evolution of angular momentum and 
grain axes subject to RATs. 

4.1. Equations of motion 

The orientation of J with respect to the magnetic field 
B is described by the following equation 



rfJ 

~dt 



— -L rad 



J 

t ■ 

L gas 



(26) 



where T ra d is the vector of radiative torque (see eq. IA1|) , 
t gas is the gas damping time. 

With the use of equation JU, we can rewrite equation 
(|26|) in the spherical coordinate system J, £, </>, assuming 
that both precession timescales of ai around J and J 
around B are smaller than the internal relaxation time: 



^-=M((H(J,£,p)) <t> )-/- 

(XL tga 



gas 



dt 



M 
1 



{F(J,£,p)U, 



(27) 
(28) 



where (H(J,£,p)) $ and (F(J, ^,p))<j> are spin-up and 
aligning components, respectively, obtained from aver- 
aging equation ([6]) and iJTJ over torque-free motion, and 
the precession angle <f> and M is given by equation ([5]). 
Also, we assume that the gas damping is isotropic, and 
we disregard the alignment by paramagnetic dissipation. 

Unlike the axisymmetric grains where the angle 9 be- 
tween ai and J is constant during torque-free motion, 
irregular grains exhibit the torque-free wobbling, and the 
conserved quantities are E, the total energy and J, the 
value of angular momentum. Therefore, for an irregular 
grain, similar to WD03, we can use p to describe the dy- 
namical evolution of grain axes. With the use of equation 
([3]) , and taking time derivative of p we obtain 



J 2 dp 
dt 



= M[h(Q u (J, £,p)>* - P J(H(J, £, 



J(p-l) fl-ph/h 



1 - h/h 



(29) 



where {QJ) is the average of torque component along the 
angular velocity given by equation ([5]), and tint is given 
by equation (TT~5|) (see Appendix B). Let us define, 

(K(J^,p))4, = (Q u (J,t,p)U/(pJ), (30) 
then, equation (|29|) becomes 



J -^ = Mp[{K{J^p)) 
(P-1) (1 



(H(J,S,p)U] 
Ph/Ii 



1 - h/h 



(31) 



Before studying grain dynamics induced by RATs, 
we first want to consider how the torque components 
{F{J,i,p))d„ (#(J,£,p))0 and (K(J,^,p))^ change with 
p and £. We adopt AMO that has the inertial property of 
triaxial ellipsoids with ratio of inertia moments h ■ h '■ 
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Fig. 6. — Torque components averaged over torque-free motion, 
(H(p)) and (K(p)), for different angles £ and for two radiation 
directions ip = 0° and 70°. These torques (H(p)) and (K(p)) coin- 
cide for p = 1, I1/I2 and I1/I3, corresponding to three stationary 
points 1, 2 and 3. 

7 3 = 2 : 1.5 : 1 and h ■ h ■ h = 1-745 : 1.610 : 0.8761. 
The later is identical with the ratio of inertia moments 
of the shape 1 (see Fig. 2]) and naturally the parameter 
q max can be changed for the AMO (LH07a). A 1.5^m 
grain size and interstellar radiation field are adopted in 
this section. 

4.2. Variation of RATs withp: stationary points for p 

In Figure [5] we plot (H(J,^,p))^ and (K(J, £,p))</> as 
functions of p for J = 107iw T , £, = 22°, 140° and tp = 
(upper panel) and 70° (lower panel) for the AMO with 
parameter q max = 1.2. 

We see that (£/"(,/, £,p))^, and (K(J,^,p)}^ decrease as 
p increases to I1/I2, and then increase as p increases to 
/1//3. This is because as p — > I1/I2, grain axes wobble 
more rigorously about J, which results in the decrease of 
torques due to averaging over wobbling with large ampli- 
tude. As p increases again, cu gets rotating stably about 
a 3 , and averaged torques increase because the wobbling 
amplitude decreases. 

In addition, {H{J,^,p)) ^ and (K(J, ^p))^ are identi- 
cal at p = 1, I1/I2 = 1-33 and I1/I3 = 2. This arises 
from the fact that at these points, oj parallel to J, so 
the projection of RATs onto J, (H(J, £,p))</,, is the same 
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Fig. 7. — Difference between (K(p)) and (H(p)) for different 
angles £ for t/> = 0° and 70° and q max = 1.2 (upper panel) and t/> = 
qo qmax _ 2 5 (lower panel). Zero points denoted by filled circles 
corresponding to (K) — (H) = are stationary points of p. In the 
upper panel, three stationary points are 1,2 and 3, corresponding 
to p = 1, I1/I2 and -Z1/-Z3. In the lower panel, in addition to 
stationary points 1,2 and 3, there exists another stationary point 
4. 



with that on u, (K(J, £,p))</>. 

The sharp decrease of (K(J,£,p)}<p when p — > I1/I2 is 
associated with the dynamics of a triaxial ellipsoid. Nev- 
ertheless, for p very close to I1/I2, w oscillates in the lab 
system, and gives rise to a sharp change of (.?r(J, £,??))</, 
(see dashed lines in Fig. [6]). 

Equation (|31|) shows that stationary points occur when 
dp/ dt = 0. In the absence of internal relaxation for large 
grains, it requires (K(J, ^,p))^ - (H(J,£,p))j, = 0. 

Figure[7|shows (K(J, £,p)}cf, — (H(J, £,p))<)i f° r the same 
AMO with q max = 1.2, V = 0°,70° (upper) and 2.5, 
ip = 0° (lower). There we omit the dependence of torques 
on J and £. We remind the reader that these cases corre- 
spond to the alignment without and with high- J attrac- 
tors, respectively, when the perfect internal alignment of 
ai with J was assumed (see LH07a). In the upper panel, 
corresponding to the alignment without high-J attrac- 
tors, there are three stationary points 1, 2 and 3. Let 
pi, P2 and P3 be the value of p at these points. It is easy 
to see that p\ = 1, P2 = I1/I2 and p$ = J1//3. In the 
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Fig. 8. — Similar to Figure[7]but for shape 1 and AMO with the 
same ratio of inertia moments Ji : I2 : I3 = 1.745 : 1.610 : 0.8761 
and q max = 2.5, and for ip = 0°. 

alignment with high- J attractors (see lower panel) , there 
exists another stationary point 4, and its value of p is 
denoted by p±. Our calculations show that P4 depends 
on q max and ip. In Figure [T2l we plot q max as functions 
of Pa for two radiation directions -0 = and 70°. 

In addition, Figure [7] shows that when p — ► h /h , the 
figure exhibits very sharp changes. In fact, we see the 
fast changes of (K (J, £,p))<f,— (H(J, £,p)A as p increases 
from slightly smaller than I1/I2 to slightly larger than 
h/h- As p = h/h, it becomes zero. 

Similar to Figure Figure [8] shows (K(J,£,p))$ — 
(H(J,£,p))<f, for the shape 1 (see FigHJ) and AMO with 
the same ratio of inertia moments I\ : h '■ h = 1.745 : 
1.610 : 0.8761 and q max = 2.5, and for ip = 0°. It can 
be seen that the AMO reproduces well the dependence 
of RATs on p for the shape 1. For instance, we observe a 
good correspondence between the stationary points ob- 
tained, but the stationary point 4 in the AMO is very 
close to the point 1 in this case (point 4 is seen in the 
lower panel of Fig. 7). 

4.3. Variation of RATs with £: stationary points for £ 

Let us consider the forms of RATs as a function of £ 
at particular values p = 1 (perfect coupling of ai with 
J) and h/h (»i -L J) for the same AMO . 

In Figure [9] we present the torque components (F(p)) 
and (H(p)} for the alignment without high-J attractor, 
i.e., q max = 1.2, V = 0° for p = 1, h/h and h/I 3 . We 
see that stationary points £ = and tt, i.e., points 1 and 
6, appear for every p. In addition to these points, forp = 
h/h and h/h, there exist two other stationary points 
at cos£ = ±9 (2,5), and cos£ = ±0.35 (3,4), respectively. 
So, in the case of wrong internal alignment, J may not 
perfectly aligned with the magnetic field. 

The upper panel of Figure [TO] is similar to that of Fig- 
ure [9l but for the case of alignment with high-J attractor 
with q max = 2.5. For p — 1, we see four stationary points 
(1,2,7,8), where (8) is a high-J attractor. Other station- 
ary points (1,3,6,8) and (1,4,5,8) are for p — h/h and 
h/h, respectively. 

4.4. Helicity versus axis of rotation 



ellipsoidal AMO, q max =1.2, ip = Q 





0.03 






~--U_J J jU J =2:1 .5:1 \ 




0.02 


; r 








0.01 


\ 


/' 


\1 


<(d) 


0.00 


= / 2 
■» 




\ 4 5 \ : 

*. .•«'»: 


V 


-0.01 
-0.02 










-0.03 




— p=',A 
p=',A 






-0.04 










1 .0 


-0.5 0.0 


0.5 1.C 



cosf 

ellipsoidal AMO, q max =1.2, f = 





0.06 








■| 'I — ?• 1 5' 1 




0.04 












0.02 










A 












CL 

_E 
V 


0.00 
-0.02 












-0.04 

-0.06 
-1 




p=1 

p=l, 

pA 


/p 






.0 


-0.5 


0.0 


0.5 1 







cosf 



Fig. 9. — Torque components (F) and (H) for p = 1, I1/I2 
and I1/I3. (1,6), (1,2,5,6) and (1,3,4,6) are stationary points for 
p = 1 , and respectively. 

For p = h/h, the rotation of the grain is about the 
axis of minor inertia a3. In this case, RATs depend on 
the angle between the grain rotation axis a3 and the ra- 
diation direction k. This angle is also the angle between 
J and B, £, when ip — 0°. Therefore, the dashed line 
in Figure [9] represents RATs as a function of the angle 
between the grain rotation axis and k. 

For p = 1, grain rotates about ai described by the an- 
gle 6. Thus, the solid line describes RATs as a function 
of = £. We can see a large difference of (F(p)) when 
the rotation axis changes from the axis of maximum mo- 
ment of inertia to axis of minor inertia. The lower panel 
is similar to the upper one, but it shows (H(p)). We 
see that the helicity of the grain for p = 1 is similar to 
that for p = h/h- I n other words, the helicity of the 
grain is identical for the rotation around the maximal 
and minimal axes. 

4.5. Trajectory map 

To visualize grain alignment, we present trajectory 
maps using J, £ and p, which are the solutions of equa- 
tions of motion ([27|) - ((3l]) . We adopt initial condition 
Jo = 4hu>T, Co spanning from to tt, and po = 1.1. Here 
we consider only one initial value po < h/h (he., grain 
initially in positive flipping state). 

In Figure [TT] we present the trajectory maps for the 
RAT alignment of the triaxial AMO induced by the ISRF 
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Fig. 10.— Similar to Fig. E] but for the case q max = 2.5, i.e. 
alignment with high-J attractor for -0 = 0°. (1,2,7,8), (1,3,6,8) and 
(1,8) are stationary points for p = 1, I1/I2 and I1/I3, respectively. 

for ip = and for q max = 1.2 and 4, corresponding to 
the alignment without and with high-J attractors (up- 
per and lower panels). In the former case, the alignment 
occurs at two low- J attractors A and B, where A cor- 
responds to the perfect alignment of J with B, and B 
corresponds to cos£ = 0.64. We see that the point B 
is the new low-J attractor appearing in the absence of 
internal relaxation. 8 . 

In the case q max = 4, the alignment has a high-J at- 
tractor A at cos£ = 1, and low-J attractors C and D 
with cos£ = —0.65 and 0.65, respectively. 

In addition, we studied the radiative alignment for the 
case q max = 0.78 and ip = 70°, which gives the trajec- 
tory map with high-J attractor. We observed that the 
alignment is similar to the lower panel in Figure 1111 

4.6. Internal alignment of grain axes with 3 

Let us discuss the internal alignment corresponding to 
the trajectory maps shown in Figure fTTl The initial value 
of p was assumed to be smaller than I\/l2- 

For the alignment without high-J attractor (the upper 
panel), the point A corresponds to the perfect internal 
alignment with p = 1, while the point B has internal 
alignment with J along a 2 axis, which corresponds to 
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8 LH07a found that for the AMO with q max 
are driven to one low-J attractor A for ip = 0° 
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Fig. 11. — Trajectory maps for the RAT alignment of a triaxial 
ellipsoid using the AMO, in the absence of internal relaxation for 
the ISRF with direction = 0° and with q max = 1.2 (upper panel) 
and q max = 4 (lower panel). A and B are low-J attractors in the 
case q max = 1.2, but there are two high-J attractors A and B for 
the case q max = 4 (lower panel; here the actual value of J for A 
and B is 10 2 times the value shown). 

P = h/h- 

For the case of alignment with high-J attractor, i.e., 
q max = 4 and ^ = qo ( lower panel of Fig rj-j^ the 

internal alignment depends on the initial value of pq. If 
Po < J34, then the lower panel of Figure [7] shows that 
(K(p)) — (H(p)) < for cos£o ~ 1- As a result, after 
some time, p decreases to p — 1, i.e., perfect internal 
alignment, and the aligning torque (F(p)) approaches 
(F(p = 1)). Therefore, the alignment returns to the case 
of perfect coupling of ai with J as studied in DW97 and 
LH07a, i.e., the alignment of J with B is determined by 
(F(p = 1)) and (H{p = 1)). LH07a found that for the 
AMO with q max > 2 and tp = 0°, the alignment has 
high-J attractors. We can see it in the lower panel of 
Figure QT] for q max = 4. Combined with the internal 
alignment, the high-J attractor corresponds to J||B and 
long grain axes perpendicular to B. For cos£o ~ — 1 ; we 
have (K(p)) — (H(p)) > 0, so p increases from p to p 4 
if po < Pa- Thus, the low-J attractor corresponds to the 
alignment with p = p^. As a result, A and B correspond 
to the internal alignment with ai||J, and C corresponds 
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Fig. 12. — Upper: Diagram showing the dependence of p4-the 
value of p at the stationary point 4 (see Fig. [7]and[8]l as function of 
q-max £ Qr £ WQ ]jg n t directions ip = and 70° . Shaded areas repre- 
sent the internal alignment when the initial value, po smaller than 
P4. Lower panel: Diagram showing the alignment type depending 
on q max and ip when po < p 4 . 

to internal alignment with ai J_ J. II po > p2, then there 
are two attractors with p = p 4 and pi- 

Since the stationary point p^ is important in determin- 
ing the type of alignment, in the upper panel of Figure 
[12] we plot q max against p 4 , for two light beam directions 
ip = 0° and 70°. Shaded areas describe the alignment 
with high-J attractor corresponding to a perfect align- 
ment of J with B and perfect internal alignment, i.e., 
ai||J. 

The lower panel of Figure [12] presents the alignment 
depending on q max and ip. 

Moreover, when the internal alignment happens with 
p = I, then, we can predict the alignment for different 
ip if we know q max . The result is presented in the lower 
panel of Figure [12] The solid curves, characterizing the 
boundary for the alignment with high-J attractor, are ob- 
tained assuming that pq < p^. We see that the alignment 
here is similar to what found in LH07a. We note that at 
p = 1, the average over the wobbling for irregular grains 
is identical with the averaging over the rotation angle (3 
used in LH07a. 

We note that the internal alignment with p — 
corresponds to the rotation of a grain along the axis of 
intermediate inertia moment a2. For torque- free motion, 



this rotation is unstable. Therefore, gas bombardment 
can destabilize this internal alignment, and grains may 
return to the alignment with short or long axes perpen- 
dicular to J. 

4.7. RAT alignment in presence of weak internal 
relaxation 

Let us consider a regime of weak internal relaxation, 
i.e., when t ra d < U nt < t gas . Figure [5] shows that grains 
with size in the range 1.5 to 6/xm correspond to this sit- 
uation. For this regime, it is convenient to introduce a 
new radiative timescale, t' rad , which is defined by (using 
eq. EH 

f > = (v£\ 

rad M\(K(p)) - (H(p))\ max > ^ ) 

where max denotes the maximal value of \(K(p)) — 
(H(p))\ as a function of p for given angles £ and ip. 

It can be seen that t' rad is about one order of magnitude 
larger than t ra d given by equation ([20| . Therefore, it 
is possible to have t' rad > tint > trad- For this case, 
one interesting effect occurs with the wrong alignment 
previously discussed is that grains are driven to right 
alignment on a characteristic time ti nt . However, for 
practical purposes, a change from t rad to t' rad is marginal 
due to a steep dependence. 

5. RAT ALIGNMENT BY DIPOLE AND QUADRUPOLE 
COMPONENTS OF RADIATION FIELD 

All earlier studies of the RAT alignment (DW96, 
DW97, WD03; LH07, LH08, HL08ab) have been done 
assuming that the radiation field can be described by a 
unidirectional beam directed along k (see Fig. [IJ with 
a degree of anisotropy 7. This assumption is consistent 
with grains near a point radiation source. In many cir- 
cumstances, e.g., in molecular clouds, and even in some 
diffuse clouds, the approximation of radiation field by 
dipole and quadrupole is more appropriate. 

Below, we study the RAT alignment induced by the 
dipole and quadrupole fields. For the sake of simplicity, 
we assume that the axis of maximum moment of iner- 
tia ai is constrained to be parallel to angular momen- 
tum J (i.e., thermal fluctuations and thermal flipping are 
completely disregarded; see DW97; LH07a). In fact, our 
study in HL08a proved that this is a sufficiently good ap- 
proximation unless we want to calculate the exact degree 
of alignment. 

5.1. Coordinate systems 

We define a lab coordinate system eie2e3 (ei-system), 
with ei parallel to B, ej 1 ei, lies in a plane formed 
by the magnetic field B and the dipole axis d for dipole 
field or symmetric axis for the quadrupole field, and e3 
perpendicular to the plane ei and e-2 (see Fig. IT51) . 9 Let 
X be the angle between d and B, and define a coordi- 
nate system e^Hde^e^ with = e3, so-called d-system. 
Now, the radiation direction k in the d-system is deter- 
mined by two angles tpk and 4>k (see Fig. 1 13|) . Since 
RATs depend on the angle between grain axes ai and 
the radiation direction k, let us define a coordinate sys- 
tem e? with ej||k, e% lying in the plane k and d, and e° 

9 The symmetric axis of the quadrupole field is defined by two 
dipole axes. 
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Fig. 13. — The lab coordinate system eie2e3 for the dipole 
field: B is the magnetic field direction, d is the dipole axis, \ is 
the angle between the dipole axis d and the magnetic field B, k is 
the radiation direction. 4>k^k describe the direction of radiation 
k in the lab coordinate system. 

perpendicular to this plane. The orientation of a grain 
in the eP system is described by three angles 0,/3 and 
$ (see the upper panel of Fig. [1] with ei are replaced 
by e?). The transformation between these coordinate 
systems are described in Appendix D. 

5.2. Dipole Component 

For a given direction of the grain in the lab coordinate 
system, the net RAT resulting from a dipole radiation 
field is given by 

Trad = 2^ Q(Wr o 2^T' ^ ' 

k=l 

where ipk is the angle between the radiation direction 
k and the dipole axis d, 7 is the degree of anisotropy, 
Ud(ipk) is the energy density of dipole field, Afifc = n/N 
is the element of solid angle in the direction ipk , and the 
summation (|33[) is performed in the range ipi- from — tt/2 
to 7r/2 (i.e., only outward radiation is accounted for). 
The spatial distribution of dipole radiation intensity is 
given by the usual expression 

Ud(ipk) = u ISRF sm 2 ip k . (34) 

We assume in equation (f3~4"|) that the energy density of 
the dipole field has the same amplitude to that of the 
ISRF, and we deal only with its spatial distribution. 

To study the RAT alignment, wc need to know 
components of torques in the lab coordinate system, 
i.e., Q ei (£,tp, <f>), Q e3 (f, V>, 4>) and Q e3 (^ip,(f>). It is 
straightforward, but tedious to obtain Q ei (£, tp, <fi) from 
<3r(Q,A $) given by equations (|A"4|) - (|Al)|) through a se- 
ries of coordinate transformation. Firstly, we perform the 
transformation from k-system to d-system, then, from 
d-system to ei-system. Finally, we average the resulting 
torques in the spherical system J, £, <f> over the azimuthal 
angles of k about dipole axis, and the Larmor pre- 
cession angle cp (see Appendix E for more detail). 

Using the obtained RATs, we solve the equations of 
motion to find J, £ as functions of time, and present their 
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Fig. 14. — Upper panel: Aligning, (F)$, and spin-up, {H)a> 
components for the AMO with q max = 3.5 induced by the dipole 
radiation field with \ = 0. Lower panel: The phase trajectory map 
obtained by torques in the upper panel. A is a high-J attractor 
corresponding to the perfect alignment of J with respect to B, and 
B is a low- J attractor. 

evolution in terms of phase trajectory maps (see LH07a 
for more detail). As an example, we use the AMO for a 
grain size a — 0.2/j,m, and q max = 3.5. For the dipole 
radiation field with \ — 0°, the torque components and 
the corresponding trajectory map are shown in the upper 
and lower panel of Figure 1141 respectively. The upper 
panel shows the existence of stationary points at cos£ ~ 
±1. It is possible to check that the stationary point 
cos£ = 1 is an attractor because (F') '(H) 1 cos^ =1 < ® 
and (H)$ > 0. This high-J attractor is denoted by a 
circle (A) in the trajectory map (see lower panel). In 
addition, there exists also a low-J attractor B, as usual. 
Both stationary points correspond to a perfect alignment 
of angular momentum with the magnetic field. 

5.3. Quadrupole Component 

Let us consider now the RAT alignment by the 
quadrupole component of radiation field. The spatial dis- 
tribution of the energy density for the quadrupole field 
is given by 

Uquad = U/SflFSi^V'fcCOS 2 ^, (35) 

where ipk is the angle between radiation direction and the 
symmetric axis of the quadrupole. Following the same 



Radiative torques alignment 



13 



procedure as with the dipole field, we obtain results for 
the AMO with q max = 0.78 (i.e., similar to the value of 
the shape 2 at A = 1.2/im, see LH07a) and the grain 
size a e ff = 0.2/zm in Figure [TBI for two directions of the 
quadrupole with the magnetic field: x = 0° (upper panel) 
and 45° (lower panel). The alignment occurs without 
high-J attractors in the former case, but with a high- J 
attractor in the later one. 

5.4. Simultaneous action of dipole and quadrupole 
components 

Figure [TBI shows the existence of high-J and low- J at- 
tractors as functions of q max and x f° r the dipole (up- 
per panel) and quadrupole (lower panel) fields. For the 
dipole field, it can be seen that the high-J attractor ap- 
pears when q max > 3.2 for x = 0°, then increases with 
X increasing, and q max > 100 for x = 30°. For x > 30°, 
the high-J attractor exists when q max < 1.8. On the 
other hand, for the quadrupole field, there is the high-J 
attractor when q max < 2.5 for a wide range x > 20°. 

LH07a have calculated RATs for a number of irregular 
grain shapes, grain size and wavelength, and found that 
for the majority of shapes we have q max < 10. Thereby, 
from Figure [TFJ we see that the possibility of alignment 
with high-J attractor is enhanced for both the dipole and 
quadrupole fields compared to a single beam. 

In the upper panel of Figure 1171 we present the ratio 
of energy density of dipole to quadrupole components, 
Ud/Ug, as a function of X for the AMO with q max = 1.2. 
For this AMO, the high-J attractor appears frequently 
(see Fig. UMivver). It reveals that the energy density of 
dipole component is required to be dominant over that 
of quadrupole in order to have high-J attractor for x < 
20°. In the range x > 30°, the high-J attractor does not 
depend on the ratio of their energy density, because both 
components produce high-J attractors. 

Apart from the existence of high- J attractors, the value 
of angular momentum J at theses points is essential for 
calculations of degree of alignment. It is easy to see that 
the angular momentum achievable (i.e., maximal value) 
depends on the angle X . 

Now, let us adopt the irregular grain (shape 1) 
with RATs obtained using DDSCAT, and calculate the 
value of the maximal angular momentum J ma x(x) = 
h^max(x) produced by dipole and quadrupole compo- 
nents of the ISRF as a function of the angle x- The 
anisotropy degree of radiation 7 = 0.1 is assumed for all 
components. We comapre the results with that induced 
by a beam of radiation. 

In the lower panel of Figure [17] we present the result- 
ing value of UJ m ax(x)- There it can be seen that uj max (x) 
decreases rapidly with increasing x for the dipole compo- 
nent, but it exhibits slow change for the quadrupole com- 
ponent. In addition, uj m ax induced by the dipole compo- 
nent is ~ 10 times greater than that by the quadrupole 
one. Therefore, when we only focus on the effect of 
spin- up by RATs, the quadrupole component can be ne- 
glected. The dipole component results in lower values of 
u m ax compared to that produced by the radiation beam 
(see Fig. [TT1 lower), as a result of the average over the 
entire space for the dipole field. 

5.5. Implications for grain alignment 
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Fig. 15. — RAT alignment by the quadrupole field with q max = 
0.78 for x = 0° (upper panel) and 45° (lower panel), respectively. 
In the upper panel, A is a high-J repellor, and B is a low- J attrac- 
tor, both corresponding to the perfect alignment of J with respect 
to B. In the lower panel, A becomes a high-J attractor and B is a 
low-J attractor. 

The degree of RAT alignment depends strongly on the 
possibility of existence of high-J attractors (see HL08a). 
In LH07a, we identified the criteria for the existence of 
high-J as a function of the angle between the beam di- 
rection and the magnetic field. There it is shown that 
the high-J attractor occurs when q max > 2 for if) < 45°, 
and q max < 1 for i/j > 45°. In other words, there is 
a gap of q max from 1 to 2, in which there is no high-J 
attractors, irrespective of grain shape and size. Fortu- 
nately, the dipole and quadrupole radiation can produce 
the high-J attractors for q max < 2 (see Fig. [TTJ) for 
X > 30°. It indicates that the dipole and quadrupole 
extends the range of high-J attractor and enable us to 
have high-J attractors for q max = 1 to 2. This range 
is easily satisfied for irregular grains (see LH07a) in a 
wide range of A/a e //. Therfore, an enhancement of the 
degree of the RAT alignment is expected for the dipole 
and quadrupole field compared to the case of a radiation 
beam. 

6. DISCUSSION 

In the rest of this section we provide a brief account of 
our accomplishments in the paper and present our out- 
look on the further work in the field of grain alignment. 
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Fig. 16. — Diagram depicting the existence of high- J and low- J 
attractors depending on q max and x f° r dipole field ( Upper panel) 
and quadrupole field (Lower panel). The shaded area represents 
the domain of the co-existence of high- J attractors and low- J 
attractors, and the white area represents the presence of low- J 
attractors only. 

6.1. Alignment in environments different from ISM: 
large grains 

Traditionally, grain alignment was the topic of inter- 
stellar medium (ISM) research. The gap between the 
studies of polarized radiation from aligned dust in en- 
vironments other then the ISM and the theory of grain 
alignment got so wide that the researchers outside the 
ISM domain sometimes write papers about aligned dust 
and do not refer to any theoretical work on grain align- 
ment. However, there is ample evidence of grain align- 
ment in non-ISM environments (see Tamura et al. 1999; 
Rosenbush et al. 2007; Hough et al. 2007). 

In many of these environments, e.g., for dust in comets 
and circumstellar dust, the radiation field is stronger 
than in the ISM, thus RATs are stronger. One difference 
that one faces there is that grains may be substantially 
larger. As we demonstrated in the paper, for larger grains 
the internal relaxation gets slower, which makes one won- 
der whether the RAT alignment happens very differently 
from the ISM case. Our study shows that the direction 
of alignment of the angular momentum J with respect 
to the magnetic field B is similar to that in the presence 
of strong internal relaxation. This alignment can be per- 
fect with high- J and/or low- J attractors, depending on 
the factor q max . However, the alignment of grain axes 
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Fig. 17. — Upper panel: Ratio of energy density of dipole to 
quadrupole components, u,i/u q , for which grains are aligned with 
high-J attractors for the AMO with q max = 1.2. The shaded area 
represents the domain of the co-existence of high-J attractors and 
low -J attractors, and the white area represents the presence of 
low- J attractors only. Lower panel: The maximal value of angular 
velocity for a grain size 0.2/im, as a function of x f° r dipole and 
quadrupole fields and one beam field. For the later, x denotes tp 
the angle between the beam direction and the magnetic field. 



with respect to J depends on the initial angle between 
them, more precisely, on the initial value of dimension- 
less parameter p. If this angle is smaller than a particular 
value, defined in the main text, which is a function of 
q max , the grain axis of maximum moment of inertia ai 
gets aligned along J, and both ai and J become aligned 
with B. HL08b showed that pinwheel torques (e.g., H2 
formation, isotropic radiative torque, etc.) increase fast 
with the grain size. As a result, large grains can be spun 
up to suprathcrmal rotation with J 3> Jth- The direct 
effect of suprathcrmal rotation is the alignment of ai and 
J. 

Table 1 compares our results for grains without and 
with internal relaxation using the AMO. We can see that 
for high-J attractors, the alignment is the same in both 
cases. We note that the results in the later case were 
shown to be consistent with the alignment of irregular 
grains obtained by DDSCAT. However the general case 
of alignment of large irregular grains in the absence of 
internal relaxation requires further studies to confirm 
the predictions obtained with the AMO. Further stud- 
ies should consider also the effect of pinwheel torques on 
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the alignment in the absence of internal relaxation. 

6.2. Utility of AMO 

The analytical model (AMO) was introduced in LH07a 
to describe the case of alignment in the presence of per- 
fect alignment of J with the axis of maximal moment of 
grain inertia. It provided excellent representation of the 
alignment of irregular grains within this approximation. 
Later in HL08a and HL08b we applied AMO to grains 
in the presence of thermal fluctuations, which at small 
values of J partially randomize the alignment of J with 
the aforementioned axis. Nevertheless, AMO happened 
to do a nice job in these cases, adequately describing the 
behavior of irregular grains. 

The case of no internal alignment is the extreme case 
of testing AMO. For instance, AMO has only one helic- 
ity axis, while multiple helicity axes are possible for an 
irregular grain. It is interesting, that even for this case, 
our limited testing is indicative of the AMO's utility. 

6.3. Towards modeling of polarized radiation 

Similar to the polarimetric work of non-ISM observers, 
the modeling of polarization arising from aligned grains 
has been mostly developing without much connection to 
the theory of grain alignment. Exceptions from this rule 
include Cho & Lazarian (2005, 2007), Pelkonen et al. 

(2007) , Bethell et al. (2007), and Falceta-Golcaves et al. 

(2008) . However, this modeling was done on the basis of 
somewhat ad hoc alignment prescriptions, which should 
be improved as the theory gets predictive. 

A step towards a more reliable polarization modeling is 
done in this paper where we considered grain alignment 
induced by the dipole and quadroupole components 1 of 
the radiation field. The earlier studies assumed that the 
radiation was coming from a particular direction, which 
is, for instance, not the case for most of the grains in 
starless cores or accretion disks. In the latter cases it is 
proper to decompose the radiation field in multipoles and 
consider the effects of the individual components. Our 
study shows that the parameter space for having high- 
J attractors differs for the alignment by the different 
multipole components. 

Moreover we identified the range of torque ratio, q max , 
for which the presence of dipole and quadroupole com- 
ponents of radiation field results in the alignment with 
high-J attractors. The required range q max is fulfilled 
for irregular grains studied in HL08a. 

6.4. Credit to Dolginov & Mytrophanov 1976 

In our paper we showed that the analytical results on 
RATs in Dolginov & Mytrophanov (1976) were incorrect 
(see Fig. [THJ). Therefore, naturally, our present results 
on the RAT alignment obtained for grains in the absence 
of internal relaxation, which was also the assumption in 
Dolginov & Mytrophanov (1976) study, differ from those 
in the study. This should not, however, undermine the 
pivotal significance of the the Dolginov & Mytrophanov 

10 The alternative to modeling the RAT alignment by the mul- 
tipoles of the radiation field is to study the alignment numerically 
at every point of the data cube using the radiation field at the 
particular point. Naturally, this would entail much more intensive 
computations. 



paper 11 . This paper is important as it discovered RATs 
and discussed the possibility of irregular grains being 
aligned by RATs, even if it failed to describe it quan- 
titatively. In addition, the notion of grain helicity, which 
is the central concept of AMO, can be traced back to 
Dolginov & Mytrophanov (1976) work. Moreover, as we 
mentioned in the introduction, the problem was so hard 
that the papers that followed Dolginov & Mytraphanov 
(1976), in all its complexity, were not able to capture 
correctly the physics of the RAT alignment either. 

6.5. Other processes 

Our paper, for the sake of simplicity, does not con- 
sider the effect of the pinwheel torques. Such torques, 
e.g. torques arising due to H2 formation over catalytic 
sites over grain surface (Purcell 1979), were discussed in 
the framework of the RAT alignment in Hoang & Lazar- 
ian (2008b) for grains with strong internal relaxation. A 
study of effect of these torques on large grains, for which 
the effects of internal relaxation are reduced, will be done 
elsewhere. In addition, in our paper we considered the 
RAT alignment in respect to magnetic field. As we dis- 
cussed in LH07a, for a sufficiently slow rate of Larmor 
rotation, the alignment can happen also in respect to the 
radiation field. If, however, grains have superparamag- 
netic inclusions the rate of rotation increases substan- 
tially. This, potentially, provides another way of testing 
whether grains have or do not have superparamagnetic 
inclusions. 

7. SUMMARY 

In this paper we continued our work on the RAT align- 
ment using both AMO and DDSCAT calculations of 
torques. Our principal results in the paper above can 
be summarized as follows: 

• We identified the range of grain size for which the 
internal relaxation within a normal paramagnetic grain 
is not important. 

• We demonstrated that, in the absence of internal re- 
laxation, RATs can align the grain's angular momentum 
with respect to the magnetic field, which is similar to the 
alignment in the presence of strong internal relaxation. 
For the internal alignment of grain axes with the angular 
momentum, it can be perfect, i.e. ai||J when the initial 
angle between them is small. 

• We studied the RAT alignment induced by dipole and 
quadrupole components of the radiation field assuming 
perfect internal alignment of grain axis of maximum mo- 
ment of inertia with the angular momentum due to strong 
internal relaxation. Using the AMO, we found that the 
parameter space for the existence of high-J attractors is 
extended compared to the earlier studied case of a sin- 
gle direction radiation. This parameter space is given 
by the range of q max and the angle x between the sym- 
metric axis of dipole and quadrupole radiation fields and 
the magnetic field. Therefore, higher degrees of radiative 
alignment are expected. 

• Our study for the joint action of dipole and 
quadrupole components showed that for the angle x be- 
tween the dipole and quadrupole axis and the radiation 

11 Incidentally, the same paper discusses for the first time the 
Barnett effect in the application to insterstellar grains. This in- 
duced the all-important notion of fast Larmor precession of grains 
and helped E. Purcell to discover the effect of Barnett relaxation. 
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TABLE 1 

Comparison of radiative alignment for AMO 



Without internal relaxation (this work) With internal relaxation (LHOTa) 

High-J attractors Low-J attractors High-J attractors Low-J attractors 

J||B J aligned parallel or J||B J aligned parallel or at 

at some angle with B at some angle with B 

Long axes _L B Long axes _L or || J Long axes _L B Long axes _L B 



RATS FOR THE ANALYTICAL MODEL: AMO 
RAT for the toy model in Figure [2] is given by 

Trad^^f^QT, (Al) 

where Qr = {Q ei ^i + Q e2 e 2 + Qe 3 ^3) is the vector of RAT efficiency, with Q ei , Q e2 and Qe 3 the components of Qr 
in the laboratory system. Here Z 2 is the size of the squared mirror, A is the wavelength, and u ra d is the energy density 
in unit erg cm' 3 of the radiation field. 

Using the self-similar scaling of the magnitude of RATs obtained for an irregular grain of size a illuminated by 
radiation field of wavelength A, 

A^ 3 



|Qr|~0.4 - for A > 1.8a, (A2) 



~ 0.4 for A < 1.8a, (A3) 
and the functional forms of RATs from the AMO, we can write RAT components as following 

o.(e,ft*-o)- '°: k rv" (6 'L* =0) . <A4) 

Qe 2 (9, 0, * = 0) = - IQrl ^^ = °\ (A5) 

0, ( e.ft.-o)- '°±7/- (e 'l^ 0) . (*») 

^J{q max ) 2 + 1 qT 3 

where 

4/i / ?7. 2 7l 2 \ 

q ei (9, 0, $ = 0) = — ±C niri 2 cos 2 9 + -^-cos/3sin29 - -^cos/3sin29 - nm 2 sin 2 9cos 2 /3 , (A7) 
A y 2 2 J 

q e2 (9, 0, $ = 0) = -^-C (n 2 cos/?cos 2 9 - ^ 2 -cos 2 /3sin29 - ^ 2 -sin29 + n^cos/3sin 2 9) , (A8) 
A V 2 2 / 

g e3 (9, /?, $ = 0) = ^-Cnisin/3 [niCos9 - n 2 cos/3sin9] + ( f] ^(s 2 - l)X(9)sin29, (A9) 
A \hj A 

with <?™ a:E = max(q ej (9, /3, $ = 0))/? for j = 1, 2 and 3. The ratio of torque components is then defined by 

n m ax = max(Q ei {e,P,$ = 0))p 

q max(Q e2 (Q,f3,<t> = 0))p' [ ) 

In equations (|A7|) - (|A9[) . C is a function defined as 

C = |ra!Cos9 — n 2 sin9cos/3| , (All) 

where 9 is the angle between the axis of maximum moment of inertia ai and the radiation direction k, (3 is the angle 
describing the rotation of the grain about ai (see Fig. [TJower); n\ — — sina, n 2 = cosa are components of the normal 
vector of the mirror tilted by an angle a in the grain coordinate system, a, 6 are minor and major semi-axes of the 
spheroid, s = a/b < 1 and e is the eccentricity of the spheroid, l\ is the distance from the mirror to the spheroid, 
and l 2 is the size of the squared mirror; K(Q) is the fitting function (see also LH07a). The second term of equation 
(|A9|) represents the torque due to the spheroid. Assuming l\ ~ A and b, a ~ I2 *C l\, then this term is subdominant 
compared to the first term (see Fig. [2]). Thus, we disregard it in our calculations. 
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FlG. 18. — Upper panei/Comparison of RATs components for AMO and DM76. is the angle between the axis of maximum moment of 
inertia (shortest) with the radiation in AMO, and is the angle between the twisted axis (longest axis) and radiation for DM76. Stationary 
points = 0,7r appear in both models, but the form of torques are different. (F) is symmetric for the AMO, but not for SM76, also in 
DM76, (H) is symmetric instead. Lower panehA brick grain adopted in calculations for internal relaxation timescale. 

Our calculations for the alignment in the presence of thermal fluctuations showed that the AMO can reproduce the 
alignment property with low-J as found with RATs obtained by DDSCAT when q ei is modified to (see HL08a) 



q ei (e,/3,^ = 0) = -^C 



nin 2 - 



[3cos 2 9 - 1] 



— cos/3sin2© -cos/3sin28 — niTi2COs2/3 



(A12) 
Also, 



This modification can arise from the imperfect scattering and/or the absorption effect by the mirror (LH07a). 
the results in LH07a remain unchanged because the averaging over (3 for the last term goes to zero. 

We adopt the AMO with a — 45° in this paper, unless mentioned otherwise. We also assume the amplitude of Q e3 
is comparable to that of Q ei and Q e2 . 

RATs at a precession angle $ (see Fig. QJower) can be derived from RATs at $ = using the coordinate system 
transformation, as follows: 



Q ei (e ) /3 ! $)=Q ei (e,A^ = o), 

Q e2 (©,/3,$) = Q e2 (e,/?,$ = 0)cos$ 
Q e3 (O, /?, *) = Q e2 (0, p, $ = 0)sin$ - 



? e3 (e,/3,$ = 0)sin$, 
? e3 (e,/3,$ = 0)cos$. 



(A13) 
(A14) 
(A15) 



AMO versus the model of Dolginov & Mytrophanov 1976 

Figure [TS] shows the comparison of the RAT components for the AMO and RATs obtained for the twisted spheroid 
from Dolginov & Mytrophanov (1976). It can be seen that the torques are radically different. According to LH07a, 
the AMO corresponds to calculations of torques for irregular grains, thus we can conclude that the model in Dolginov 
& Mytrophanov (1976) does not represent adequately RATs. 

DYNAMICAL EQUATIONS FOR IRREGULAR GRAINS 
For irregular grains, we define a parameter which is constant during the free-torque motion: 

2hE 



V 



J 2 



(Bl) 



where E is the total energy, I\ is the inertia moment along the principal axis ai and J is the value of angular momentum. 
Since both J and E are conserved during the torque-free motion, p is accordingly conserved. The evolution of p in 
time due to external torques is then 



dp 
~dl 



The energy of the grains varies as 



1 

= J2 

dE 
~dT 



J ^ 
'~dt 



2I X E 



T.u 



dJ 



(B2) 



(B3) 
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Fig. 19. — Upper panel: RAT alignment by Q e ±, in the absence of internal relaxation. A and D denoted by circles are high- J attractors, 
and C and D denoted by crosses are repellors. The alignment at A corresponds to ai X B, and ai||B for the alignment at D. Lower panel: 

SThfiuSs^ iPBr&id^ nuclear relaxation. The second terms represents 

the effect of external torques on the grain rotational energy. This equation differs from that in WD03 by a factor of 
2. We consider here only RATs, so T = T ra d — MQy. 
Substituting dE/dt into equation (|B2jl . we get 



dp 1 r ^ / - du> 
dt = jz^ 

Averaging equation (|B4[) over the torque-free motion, we obtain 



h j.- + r..)-pj.r 



(B4) 

(B5) 
(B6) 



dt u dt t^t V 1 - h/h 

where we used T = dJ/dt, and 

Qui = Qr-^ = Qoi^l + Qa 2 U2 + Qa 3 U 3 - 

Here u>i for i=l,2 and 3 corresponds to the components of angular velocity along three grain axes ai,a2 and a3. 
RAT ALIGNMENT INDUCED BY A SINGLE COMPONENT FOR THE SPHEROIDAL AMO 

We study first the alignment with high-J attractor. For this purpose, we consider the alignment by the first 
component Q ei , which is shown to produce the alignment with high-J attractor (see LH07a). With this simplification, 
we can obtain the analytical results for the motion. According to equation (31) in HL08a, the averaged value of Q ei 
over torque-free motion, and fast Larmor precession is given by 

(Qe^^QTr ( 3cosWc?+ 3sin2 f n ^ - , (CI) 



(Q ai )^ = Q™ a:E cos6>cosC ( -1 + 3cos 2 £cos 2 (9 



9 • 2t ■ 2, 

-sin fsin t 
2 



(C2) 



The aligning and spin-up torque components are then 

CF% = -(Q ei Usin£, (C3) 
(H)* = (Q ei )0cose (C4) 

It can be seen that the torque components are functions of two alignment angles 9 and £. 

Let us investigate the property of these stationary points. Consider first the stationary points £ = ir and 9 = 90°. 
Equation (JCU) shows that (iJ) = (Q™ ax ) 2 > because (Q ei ) = -Q™ ax . In addition, the first derivative of (F)^ < 
(see eq. IC3|) . Therefore, £ = tt is a high-J attractor with 9 = 90°. Moreover, the stationary point £ = 0, 9 = is still 
a high-J attractor because {F')^/ (H)$ < 0. This indicates that in the absence of internal dissipation, the alignment 
occurs in two types: the longest axis parallel and perpendicular to the magnetic field. The former is consistent with 
the Davis-Greenstein prediction, i.e., "right" alignment, while the later provides "wrong" alignment (see LH07a). 

The trajectory map in Figure fT9l for the alignment by one component Q e i exhibits two high-J attractors A and D, 
as expected. In addition, there are two repellors B and C. The final state of grains depends on their initial angles 
£o- For instance, grains with initial angles |cos£o| > 0.6 are aligned on A and D, but grains with |cos£o| < 0.6 are 
constrained within two repellors B and C, and finally damped by gas friction. 

The lower panel represents the evolution of 9 and £ as functions of time for few grains with cos£o < —0.8. It can 
be seen that the angular momentum of these grains are perfectly aligned with respect to the magnetic field, but the 
shortest axis ai is nearly perpendicular to J, corresponding to "wrong" internal alignment. 
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TRANSFORMATION OF COORDINATE SYSTEM 

To find the torques in the lab system (e;) when the radiation direction k varies, we need to implement the coordinate 
transformation from the system kk x k y to ei (see Fig. [13]). We assume that the magnetic field is directed along ei, 
then the kj = k, and k y lies in the plane B, k. Thus, k x is perpendicular to the plane k, k y . Therefore, the coordinate 
system k, k x , k y acts as the system ei in Figure [TJ 

Denote Qk be torques components in the k-system, we need to know the torque components in the lab system. The 
transformation from k-system to ej-system is carried out by three rotations with three Euler angles. 

Following Figure I13[ we have 

= eicosx + e 2 sinx, (Dl) 
e\ = — eisinx + e2COSX, (D2) 
e* = e 3 (D3) 

The radiation direction k is determined by two angles ?pk and <pk hi the lab system, we have the coordinate trans- 
formation from k-system to d-system: 

e i = eJcosV'k + e^sin^kcos^k + egsin^ksirx^k (D4) 

e% = — e^sin-0k + e^cos-^kcos^k + egCOS-fAkSin^k (D5) 

e 3 = e^sin^k + e^cos^ (D6) 
Assuming that ai is parallel to J, which is described by £, <f> in the lab system, then we have 

ai = eicos^ + e2sin£cos0 + e3sin£siri0 (D7) 
a 2 = — eisin^ + e 2 cos£cos^> + e 3 cos£sin</>, (D8) 

(D9) 

We can find the angles 0, $ and /3 in the e? system: 

cos9 = e?.a x , (D10) 

$ sin6-ai.e2 

tan- = Dll) 

2 ai.eg v ; 

i / sinG + a^.eS \ 

P = 2tan —7T( o ^ o • ^ . D12 

\smfc)(a2.e3Cos<P — a2.e2Sin$) / 

By substituting equations ([D4 jl -([D6 ]l and (jD7| and (|D8| into (jDlOj) and (|D12|) . we obtain the angles 9, $ and (3, the 
torques Q° ei are interpolated. Finally, we find the corresponding torques, Qei(£,,(j>) by coordinate transformation: 

Qe t =Cf 3 Bj k Ql k , (D13) 

Here matrices C T ,B T are tranposal matrices of C and B, inferred from the equations (|D10|) - (|D12[) and (|A4|) - (|F31|) 
following 

ef = Cyej, (D14) 
e? = Bye* (D15) 



